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Abstract 

We apply a multi-color extension of the Beck-Fiala theorem to show that the multiobjective 
maximum traveling salesman problem is randomized 1/2-approximable on directed graphs and 
randomized 2/3-approximable on undirected graphs. Using the same technique we show that the 
multiobjective maximum satisfiablilty problem is 1/2-approximable. 



1 Introduction 
> 

en . We study multiobjective variants of the traveling salesman problem and the satisfiability problem. 

[~^ . • The ^-objective maximum traveling salesman problem: Given is a directed/undirected com- 

^^ I plete graph vi^ith edge weights from N'^. Find a Hamiltonian cycle of maximum weight. 

• The ^-objective maximum weighted satisfiability problem: Given is a Boolean formula in 
conjunctive normal form and for each clause a non-negative weight in N . Find a truth 
^ . assignment that maximizes the sum of the weights of all satisfied clauses. 



In general we cannot expect to find a single solution that is optimal with respect to all objectives. 
Instead we are interested in the Pareto set which consists of all optimal solutions in the sense that 
there is no solution that is at least as good in all objectives and better in at least one objective. 
Typically, the Pareto set has exponential size, and this particularly holds for the traveling sales- 
man and the satisfiability problems considered here. We are hence interested in computing an 
approximation of the Pareto set. 

A popular strategy for approximating single-objective traveling salesman and single-objective sat- 
isfiability is to compute two or more alternatives out of which one chooses the best one: 

• For each cycle in a maximum cycle cover of a graph, remove the edge with the lowest weight, 
and connect the remaining paths to a Hamiltonian cycle. 
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• For some formula, take an arbitrary truth assignment and its complementary truth assign- 
ment, and return the one with the highest weight of satisfied clauses. 

However, in the presence of multiple objectives, these alternatives can be incomparable and hence 
we need an argument that allows to appropriately combine incomparable alternatives. 

While previous work focused on problem-specific properties to construct solutions of good quality, 
we show that the Beck-Fiala theorem [BF81J from discrepancy theory and its multi-color extension 
due to Doerr and Srivastav |DS03j provide a general and simple way to combine alternatives 
appropriately. Its application leads to simplified and improved approximation algorithms for the 
fc-objective maximum traveling salesman problem on directed and undirected graphs and the k- 
objective maximum weighted satisfiability problem. 



2 Preliminaries 

Multiobjective Optimization Let k > 1 and consider some ^-objective maximization problem 
O that consists of a set of instances I, a set of solutions S{x) for each instance x € Z, and a 
function w assigning a fc-dimensional weight w{x, s) € N to each solution s G S{x) depending also 
on the instance x G Z. If the instance x is clear from the context, we also write w{s) = w{x,s). 
The components of w are written as Wi. For weights a = (ai, . . . , a^), b = (61, . . . , 6^) € N^ we 
write a >b if Ui > bi for all i. 

Let X S X. The Pareto set of x, the set of all optimal solutions, is the set {s S S{x) \ -<3s' G 
S{x) (w;(x,s') > ■w{x,s) and w{x,s') ^ w(x,s))}. For solutions s,s' £ S{x) and q < 1 we say s 
is a-approximated by s' if Wi{s') > a ■ Wi{s) for all i. We call a set of solutions a-approximate 
Pareto set of x if every solution s £ S{x) (or equivalently, every solution from the Pareto set) is 
a-approximated by some s' contained in the set. 

We say that some algorithm is an a-approximation algorithm for O if it runs in polynomial time 
and returns an a-approximate Pareto set of x for all inputs x S Z. We call it randomized if it 
is allowed to fail with probability at most 1/2. An algorithm is an FPTAS (fully polynomial-time 
approximation scheme) for O, if on input x and < e < 1 it computes a (1 — e)-approximate Pareto 
set of X in time polynomial in i/e + length(x). If the algorithm is randomized it is called FPRAS 
(fully polynomial-time randomized approximation scheme). If for each e, the runtime of such an 
algorithm is polynomial in length(x), we call it PTAS (polynomial-time approximation scheme) or 
PRAS (polynomial-time randomized approximation scheme). 

Graph Prerequisites An ^''-labeled directed (undirected) graph is a tuple G = (V,E,w), where 
V is some finite set of vertices, E C V x V [E Q (2)) is a set of directed (undirected) edges, and 
w: E ^ 'N^ is a /c-dimensional weight function. If E = {V x V) \ {{i,i) \ i £ V} {E = (^)) then 
G is called complete. We denote the i-th component of w by Wi and extend w to sets of edges 
by taking the sum over the weights of all edges in the set. A cycle (of length ?n, > 1^ in G is an 
alternating sequence of vertices and edges t;o, ei, vi, . . . e^, 'Vmi where Vi G V, Cj G E, Cj = {vj-i,Vj) 
(cj = {vj-i,Vj}) for all < i < m and 1 < j < m, neither the sequence of vertices vo,vi, . . . , Vm-i 



nor the sequence of edges ei, . . . ,em contains any repetition, and Vm = vq. A cycle in G is called 
Hamiltonian if it visits every vertex in G. A set of cycles in G is called cycle cover if for every 
vertex f € 1^ it contains exactly one cycle that visits v. For simplicity we interpret cycles and cycle 
covers as sets of edges and can thus (using the above mentioned extension of w to sets of edges) 
write w{C) for the (multidimensional) weight of a cycle cover C of G. 



Approximating Cycle Covers We will consider approximation algorithms for the multiobjec- 
tive traveling salesman problem using a multiobjective version of the maximum cycle cover problem. 
For directed input graphs we have the following problem definition. 

fc-Objective Maximum Directed Edge-Fixed c-Cycle Cover (fc-c-MaxDCCp) 

Instance: N'^-labeled complete directed graph (V, E, w) and F C E 
Solution: Cycle cover C (^ E with at least c edges per cycle and F Q G 
Weight: w{C) 

For undirected input graphs we analogously define the ^-objective maximum undirected edge- 
fixed c-cycle cover problem (fc-c-MaxUCCF, for short). Let fc-c-UCC (fc-c-DCC) denote the 
problems we obtain from fe-c-MaxDCCp (fe-c-MaxUCCF) if we require F = 0. Using this notation 
we obtain the usual cycle cover problems fc-MaxDCC as /c-O-DCC and fc-MaxUCC as A:-0-UCC. 

Manthey and Ram |MR09j show by a reduction to matching that there is an FPRAS for A;-objective 
minimum cycle cover problems. The same technique can be used to show that there are FPRAS 
for fe-MaxDCC and fc-MaxUCC |Man0 9]. We show that there are FPRAS for fc-2-MaxDCCF and 
fe-S-MaxUCCp by a reduction to A;-MaxDCC and A;-MaxUCC. 

Theorem 2.1. For every k>l, A;-2-MaxDCCF and A;-3-MaxUCCF admit an FPRAS. 

Proof. For every / > 1, let /-MaxDCC-Approx (/-MaxUCC-Approx) denote the FPRAS for l- 
MaxDCC (Z-MaxUCC). We begin with the directed case. 

Let A; > 1. On input of the N'^-labeled complete directed graph G = {V,E,w) and F C E, let 
G' = (y, E, w'), where w' : E -^ N''+^ such that for all eG E, 

'w^{e) = Wi{e) for 1 < i < k 

, Jl ifeGF 

w^,i[e) = < 

I otherwise. 

For e > 0, apply {k + l)-MaxDCC-Approx to G' with approximation ratio e' = min{e, "^/{r+i)}, 
where r := ^F and return the obtained set of cycle covers that contain all edges from F. 

Let G be some (arbitrary) cycle cover with F Q G. If no such cycle cover exists, we are done. 
Otherwise, we have w'f^,^{G) = r, and with probability at least 1/2 the FPRAS must have returned 
some cycle cover C" that e'-approximates G. By e' < y{r+i) we have u^^_,_i(C") > (1 — e') ■ 
'W'fc+i(C') > (1 — V(''+i)) ■ r = r — "^/(r+i) > r — 1 and hence F C C". Moreover, by e' < e we have 



Wi{C') = w'iiC) > (1 - e') ■ w[{C) >(!-£)• w'iiC) = {I - e) ■ Wi{C) for all 1 < i < A;. Since an 
arbitrary cycle in a complete directed graph has length at least two, the assertion is proved. 

The proof for the undirected case is very similar, as we call (/i: + l)-MaxUCC-Approx instead. Since 
in a complete undirected graph every cycle has length at least three, the assertion follows. D 



Boolean Formulas We consider formulas over a finite set of propositional variables V ^ where 
a literal is a propositional variable v ^ V or its negation v, a clause is a finite, nonempty set of 
literals, and a formula in conjunctive normal form (CNF, for short) is a finite set of clauses. A 
truth assignment is a mapping I : V ^ {0, 1}. For some v (zV, we say that / satisfies the literal v 
if I{v) = 1, and / satisfies the literal v if I{v) = 0. We further say that / satisfies the clause C and 
write /(C) = 1 if there is some literal I € C that is satisfied by I. 



3 Multi-Color Discrepancy 

Suppose we have a list of items with (single-objective) weights and want to find a subset of these 
items with about half of the total weight. The exact version of this problem is of course the NP- 
complete problem partition |GJ79j . and hence it is unlikely that an exact solution can be found 
in polynomial time. If we allow a deviation in the order of the largest weight, this problem can 
be solved in polynomial time, though. Surprisingly, this is still true if the weights are not single 
numbers but vectors of numbers, which follows from a classical result in discrepancy theory known 
as the Beck-Fiala theorem |BF81j . It is important to note that the allowed deviation is independent 
of the number of vectors since this enables us to use this result in multiobjective approximation 
for balancing out multiple objectives at the same time with an error that does not depend on the 
input size. 

In the Beck-Fiala theorem and the task discussed above, we have to decide for each item to either 
include it or not. In some situations in multiobjective optimization, though, a more general problem 
needs to be solved: There is a constant number of weight vectors for each item, out of which we have 
to choose exactly one. Doerr and Srivastav |DS03j showed that the Beck-Fiala theorem generalizes 
to this so-called multi-color setting with almost the same deviation. Their proof implicitly shows 
that this choice can be computed in polynomial time. For completeness we restate the proof and 
argue for polynomial-time computability. 

For a vector x € Q™ let ||x||oo = maxj \xi\, and for a matrix A G Q^^i^ let ||^||i = maxj ^^ Wijl- 
For c > 2, n > 1 let M^ = {x G (Q n [0,1])^=" | Efe=o^c6-fc = 1 for all b € {!,... ,n}} and 
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Mc,n = Me,„ n {0, 1}^ 

Theorem 3.1 (Doerr, Srivastav [D S03] ). There is a polynomial-time algorithm that on input of 
some A G Q™xcn^ m,n G N, c > 2 andp G Mc^n finds a coloring % £ ^c,n such that \\A{p — x) I loo < 
2p||i. 

Proof. Let A := ||^||i. We start with x = X = p G Mc^n and will successively change it to a 
vector in Men- We will first describe the algorithm and then argue about its runtime. 



Let J := J(x) := {j G {l,...,cn} [ Xj ^ {0,1}} and call the columns from J floating. Let 
/ := I{x) '■= {i G {l,...,m} \ J2jeJM\^v\ ^ 2A}. We will ensure that during the rounding 
process the following conditions are fulfilled (this is clear from the start, because x = p)'- 



{A{p-x))\i = (CI) X^Mcn (C2) 

Let us assume that the rounding process is at step t where the current coloring is % = x and the 
conditions (CI) and (C2) hold. If there is no floating column, i.e., J = 0, then x ^ -^c,n and thus 
X has the desired form. 

Otherwise, assume that there are still floating columns. Let i? = {6 € {1, . . . , n} | 3A; € {0, . . . , c — 
1}: cb — k £ J} he the c-blocks that contain floating columns. Since x ^ Mc,m a c-block of x 
contains either none or at least two floating columns, thus ^B < 2^^^- 

Since 

m 

#J • A = j; A > j; j; |a.,| > 2^2^ |a.,| = EE l^^^'l > E ^^ = #^ " ^^ 
jeJ j&J i=i j&J i£i iei jeJ iei 

it holds that #/ < 2^'^- Consider the inhomogeneous system of linear equations 

{A{p-x))\i = 

c-l 

^ Xcb-k = 1 for 6 G B 



A;=0 



where each Xj is considered as a variable if j € J and as a constant if j ^ J. This system consists of 
at most #1 + #B < ^#J + ^#J = #J equations and #J variables and hence is under-determined. 
Note that the system has the solution xij because x fulfills the conditions (CI) and (C2). Since it 
is under-determined, it also has a second solution x G Q'^ . We extend x to xe G Q'^" by 



ixE)j 



Xj if j € J 
Xj otherwise. 



Consider the line {(1 — A) x + Ax^; | A G Q}. Each point on this line (or rather its restriction 
to the components in J) fulfills the system of equations and thus condition (CI). By condition 
(C2) and the definition of J it holds that < Xj < 1 for all j (z J and thus there is some A G Q 
such that X '■= (1 — A)x + Ax^ G Mc^n and at least one component becomes or 1, i.e., 

j(^(t+i)) c J(xW). Note that x^*"^^^ fulfills (CI) and (C2) even for the larger sets J(x^*^) and 
I{x^^')- Continue the rounding process with x •= X • 

Since at least one column is removed from J in each iteration, the rounding process will eventually 
stop. Let X be the final value of the coloring. We show ||yl(p — x)||oo < 2A. Let \ < i < m. 
Since at the end, J = we also have / = 0. Let x b^ the first coloring such that i ^ I. Since 
X^*^ fulfills (CI) also for /(x^*'"^^) (or /(x*^°^) if t = 0) we have {A{p - x^*)))* = 0. Furthermore it 



holds that x)'' = Xj ^ov all j ^ J{x ) ^iid \x- ~ Xj\ < 1 foi" ^-H i ^ J{x )• Finally note that 






E 



■ieJ(x(*)) l"'i 



< 2 A since i ^ I{x )• Combining these facts, we obtain 



Uip - xm = \{A{p - x^'% + {Aix^'^ - x)\ 



ieJ(x(*)) 



(t) 



Xi)l<2A. 



We now analyse the runtime. Note that we have at most en iterations, which is polynomial in the 
input length. In each iteration, we have to solve an inhomogeneous system of linear equations and 
we have to find a certain A € Q. The system, whose size is polynomial in the input length, can be 
solved in polynomial time (see for instance |GLS88^ Theorem 1.4.8]). By adding an equation of the 
form Xj = 2 for some suitable j G J, we can find a solution different to x- The value for A can be 
obtained in polynomial time by successively trying to fix each floating column to or 1, solving for 



A and checking if the resulting vector is still in M^^n- 



a 



Corollary 3.2. There is a polynomial-time algorithm that on input of a set of vectors v^'^ G Q™' for 
l<j<^, l<r<c computes a coloring x '■ {^, . . . , n} — )• {1, . . . , c} such that for each 1 < i < m 
it holds that 



^ n c n 



j.x(i) 



=1 r=l 



< 2?7imax \vl'^\ 

j,r 



Proof. The result is obvious for c = 1. For c > 2, we use Theorem 13. 11 Because the error bound is 
different for each row, we need to scale the rows of the vectors. Let 5i = max 



■J,r 



M'^l for 1 < i < m. 



Let A 



"1,3 



e 



where a. 



'«>(c(i-i)+'') 



hr (if s^ 



0, set it to 0) and p £ 



such 



that Pi = - for all 1 < i < en. We obtain a coloring x ^ {0, 1}^" such that for each 1 < j < n 
there is exactly one 1 < r < c such that Xc{j-i)+r = 1 ^^'^ if holds that \\A{p — x)\\oo < 2||A||i. 
Note that because of the scaling, the largest entry in ^ is 1 and thus we have ||A||i < m. Define 



X' : {1, . . . , ra} -^ {1, . . . , c} by x'ij) 



Xc{j — l)+r 



1. We obtain for each 1 < i < m: 



2m5i>2\\SiA\\i > {{Si Aip - x))i 



"^diaij'ipji -Xj' 
i'=i 



EE>r-E-^ 



j.x'(j) 



j=l r=l 



i=i 



D 



4 Approximation of Multiobjective Maximum Traveling Salesman 



Definition Given some complete N'^-labeled graph as input, our goal is to find a Hamiltonian cycle 
of maximum weight. For directed graphs this problem is called /c-objective maximum asymmetric 
traveling salesman (A;-MaxATSP), while for undirected graphs it is called /c-objective maximum 
symmetric traveling salesman (/c-MaxSTSP). Below we give the formal definition of /c-MaxATSP, 
the problem A;-MaxSTSP is defined analogously. 



fe-Objective Maximum Asymmetric Traveling Salesman (fc-MaxATSP) 

Instance: N'^-labeled directed complete graph {V, E, w) 
Solution: Hamiltonian cycle C 
Weight: w{C) 

Previous Work In 1979, Fisher, Nenihauser and Wolsey |FNW79] gave a Y2-approxiniation algo- 
rithm for the single-objective maximum asymmetric traveling salesman (1-MaxATSP) by removing 
the lightest edge from each cycle of a maximum cycle cover and connecting the remaining paths to 
a Hamiltonian cycle. Since undirected cycles always contain at least three edges, this also showed 
that the single-objective maximum symmetric traveling salesman (1-MaxSTSP) is 2/3-approximable. 
Since then, many improvements were achieved, and currently, the best known approximation ratios 
of 2/3 for 1-MaxATSP and 7/9 for 1-MaxSTSP are due to Kaplan et al. |KLSSn5j and Paluch, Mucha 
and Madry jPMMn9j . 

Most single-objective approximation algorithms do not directly translate to the case of multiple 
objectives, and hence we need more sophisticated algorithms. For A;-MaxATSP and fc-MaxSTSP, 
where k > 2, the currently best known approximation algorithms are due to Manthey, who showed 
a randomized (1/2 — e)-approximation of /c-MaxATSP and a randomized (2/3 — e)-approximation of 
/c-MaxSTSP |Man09| . Recently, Manthey also showed a deterministic (V2fe — e)-approximation of 
A;-MaxSTSP and a deterministic (Y(4A;-2) — e)-approximation of /c-MaxATSP [Manllj . 



Our Results We show that A;-MaxATSP is randomized 1/2-approximable and A;-MaxSTSP is 
randomized 2/3-approximable using the following idea. We choose a suitable number / depending 
only on k and try all sets of at most I edges F using brute force. For each such F we apply the 
FPRAS for A;-2-MaxDCCF (A;-3-MaxUCCF), which exists by Theorem EH fixing the edges in F. 
For all cycle covers thus obtained, we select two (three) edges from each cycle and compute a 
2-coloring (3-coloring) of the cycles with low discrepancy with regard to the weight vectors of the 
selected edges. Using this coloring, we remove exactly one edge from each cycle and connect the 
remaining simple paths to a single cycle in an arbitrary way. Since the coloring has low discrepancy, 
we only remove about one half (one third) of the weight in each objective. The introduced error is 
absorbed by choosing suitable heavy edges F at the beginning. The described procedure generally 
works for arbitrary c-cycle covers. 

Lemma 4.1. Let c > 2 and k > 1. If there exists an FPRAS for fc-c-MaxDCCp CA:-c-MaxUCCF, 
resp.), then the algorithm Alg-A:-MaxTSP computes a randomized (1 — ^/c) -approximation for k- 
MaxATSP ("yfc-MaxSTSP, resp.). 

Proof. Let k > 1, c > 2, and G = {V,E,'w) be some N -labeled (directed or undirected) input 
graph with m = ^V sufficiently large. 

We will first argue that the algorithm terminates in time polynomial in the length of G. Since 
there are only polynomially many subsets Fh,Fl C E with cardinality bounded by a constant, 
the loop in line [J is executed polynomially often. In each iteration the FPRAS on G = {V, E, w) 
and Fh U Fl (^ E terminates in time polynomial in the length of G and Fh U Fe, which means 
that the set S contains only polynomially many cycle covers. Hence, for each iteration of the loop 
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Algorithm: Alg-k-MaxTSP (V, E , w) with parameter c > 2 



Input : N^-labeled directed/undirected complete graph G = {V, E, w) 
Output: set of Hamiltonian cycles of G 

1 foreach Fh,Fl C E with #Fh <3ck^, #Fl < c#Fh do 

2 let 5 S N'^ with 6i = max{n S N | there are 3ck edges e € Fh with Wi{e) > n}; 

3 foreach e (z E\ Fh do 

4 if w{e) ^ 5 then modify w such that w{e) = for the current iteration of line[Tl 

5 compute (1 - V#V')-approximation S of fc-c-MaxDCCF / fc-c-MaxUCCF on (G, FhUFl) ; 

6 foreach cycle cover S" G 5 do 

7 let Ci, . . . , C^ denote the cycles in S; 

8 if for each i G {1, . . . , r}, Cj \ -F/f contains a path of length c then 

9 foreach i G {1, . . . , r} do choose path ej^i, . . . , ej^c G Cj \ F/^ arbitrarily; 
10 compute some coloring x'- {1; • • • ; '^l ^ {1) • • • ; c} such that 



r ^ r c 



c 



and remove the edges {ej^;^(j) | i = 1, . . . , r} from 5; 
11 output the remaining edges, arbitrarily connected to a Hamiltonian cycle; 



in line [H the loop in line [6] is also executed at most polynomially many times, and overall we 
have polynomially many nested iterations. In each nested iteration where each cycle of the cycle 
cover contains a path as required, we compute a coloring of {1, ■ ■ ■ ,r} with low discrepancy. By 
Corollary 13.21 this can be done in polynomial time. Observe that all further steps require at most 
polynomial time, and hence the algorithm terminates after polynomially many steps. 

Next we argue that the algorithm will succeed with probability at least 1/2. Observe that the 
only randomized parts of the algorithm are the calls to the randomized cycle cover approximation 
algorithm in line El Using amplification we can assume that the probability that all the calls to 
this algorithm succeed is at least 1/2. 

It remains to show that if the algorithm Alg-A;-MaxTSP succeeds, it outputs some (1 — i/c)- 
approximate set of Hamiltonian cycles. Hence, for the remainder of the proof, let us assume 
that the algorithm and hence all calls to the internal FPRAS succeed. Furthermore, let R C E 
be some Hamiltonian cycle of G. We will argue that there is some iteration where the algorithm 
outputs an (1 — i/c)-approximation of R. 

For each 1 < i < k, let Fn^i C i? be some set of Sc/c heaviest edges of R in the i-th component, 
breaking ties arbitrarily. Let Fh = Uj=i PH,i- We define F^ C R such that F^ n Fh = and each 
edge in Fh is part of a path in F^ U Fh that contains c edges from F^. This is always possible as 
long as R is large enough. We now have #-Fff < 3ck^ and i^F^ < c^Fh- Hence in line [J there 
will be some iteration that chooses Fh and Fl- We fix this iteration for the remainder of the proof. 

Let (5 G N as defined in line [2] and observe that 6i = mm{wi{e) \ e G Fn^i} for all i, which means 
that for all edges e G R\Fh we have w{e) < 5. Hence the loop in line [3] sets the weights of all edges 



e £ E\R that do not fulfill w{e) < S to zero, and these are the only weights that are modified. In 
particular, this does not affect edges in R, hence w^R) remains unchanged. Note that since we do 
not increase the weight of any edge and do not change the weight of the edges in R, it suffices to 
show that the algorithm computes an approximation with respect to the changed weights. 

Next we obtain a (1 — Y#v)-approximate set S of c-cycle covers of G that contain FhUFl. Since R 
is a c-cycle cover of G with FhUFl C i?, there must be some c-cycle cover S £ S with FhUFl C S 
that (1 — i/#v)-approximates R. Hence in line [6] there will be some iteration that chooses this S. 
Again we fix this iteration for the remainder of the proof. 

As in line El let Gi, . . . ,Cr denote the cycles in S. Note that each cycle contains at least c edges. 
Since each edge in Fh is part of a path in Fh U Fl with at least c edges from Fl, we even know 
that each cycle contains at least c edges not from Fh and thus the condition in line [8] is fulfilled. 
Let these edges eij be defined as in the algorithm. Note that since eij ^ Fh we have w{eij) < 5 
for all i,j, because the weight function was changed accordingly. 

In line [lO] we compute some x ^ {1; • • • ) ''I ~^ {l; • • • ; c} such that 

r ^ r c 

i=l i=l jf=l 

< 2k-6 + --w(S\FH). 
c 

Recall that by Corollary 13.21 such a coloring exists and can be computed in polynomial time. 
Removing the chosen edges breaks the cycles into simple paths, which can be arbitrarily connected 
to a Hamiltonian cycle R'. For the following estimation note that 6 < "'i ^-^ and w{Fh) > ^^w{R) 
and recall that m = #V = #R. 

r 

w{R') > w{S) -^w{ei^^^i^) 

> w(S)-2k-6---w(S\FH) 

c 

--]w(S) + -w(FH)-2k-6 
cj c 

> fi-l.\^^s) + ^w{FH) 
\ cJ 6c 

> f 1 - 1 Vl - 1^ wiR) + -w{R) 
\ cJ \ mj m 

i-i).(fl)+(-(i-l)+t)i»(i;) 

> |l-i)a.(R) 



C^ 

This proves the assertion. D 

It is known that 1-c-DCC is APX-hard for all c > 3 |BMn5j and that 1-c-UCC is APX-hard for 
c > 5 |Man08] . This means that, unless P = NP, there is no PTAS for these problems |CK99] (and 



especially not for the variants with fixed edges). Furthermore, the existence of an FPRAS or PRAS 
for these problems implies NP = RP and thus a collapse of the polynomial-time hierarchy, which 
is seen as follows. 

If an APX-hard problem has a PRAS, then all problems in APX have a PRAS and hence MAX- 
3SAT has one. There exists an e > and a polynomial-time computable / mapping CNF formulas 
to 3-CNF formulas such that if x € SAT, then f{x) G 3SAT; and if x ^ SAT, then there is no 
assignment satisfying more than a fraction of 1 — e of /(x)'s clauses |AL971 Theorem 10.1]. The 
PRAS for MAX-3SAT allows us to compute probabilistically a (1 — '^/2)-approximation for f{x) 
which in turn tells us whether or not x G SAT. Since this procedure has no false negatives we get 
RP = NP, which implies a collapse of the polynomial-time hierarchy jLau831 |Sip83| . 

So it seems unlikely that there is a PRAS for 1-c-DCC where c > 3 and 1-c-UCC where c > 5. 
However, this does not necessarily mean that the above algorithm is useless for parameters c > 3 
in the directed and c > 5 in the undirected case: The algorithm could still benefit from a constant- 
factor approximation for A;-c-MaxUCCF or fc-c-MaxDCCp. A simple change in the estimation shows 
that if the cycle cover algorithm has an approximation ratio of a, the above algorithm provides an 
approximation with ratio a(l — i/c). 

Theorem 4.2. Let k>l. 

1. /c-MaxATSP is randomized ^/2-approximahle. 

2. A:-MaxSTSP is randomized "^/s-approximable. 

Proof. We combine Theorem 12.11 and Lemma 14. li D 

5 Approximation of Multiobjective Maximum Satisfiability 

Definition Given a formula in CNF and a function that maps each clause to a /c-objective weight, 
our goal is to find truth assignments that maximize the sum of the weights of all satisfied clauses. 
The formal definition is as follows. 

fc-Objective Maximum Weighted Satisfiability (fc-MaxSAT) 

Instance: Formula H in CNF over a set of variables V, weight function w: H ^ N 

Solution: Truth assignment I: 1/ — )• {0, 1} 

Weight: Sum of the weights of all clauses satisfied by I, i.e., 'w{I) = Yl ^(C*) 

CeH 
I{C)=l 

Previous Work The first approximation algorithm for maximum satisfiability is due to John- 
son [Joh74| . whose greedy algorithm showed that the single-objective 1-MaxSAT problem is 1/2- 
approximable. Further improvements on the approximation ratio followed, and the currently best 
known approximation ratio of 0.7846 for 1-MaxSAT is due to Asano and Williamson |AW02j . 
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Only little is known about k-MaxSAT for k > 2. Santana et. al. [SBLLOO] apply genetic algorithms 
to a version of the problem that is equivalent to /c-MaxSAT with polynomially bounded weights. 
To our knowledge, the approximability of /c-MaxSAT for k > 2 has not been investigated so far. 

Our Results We show that /c-MaxSAT is 1/2-approximable mainly by transferring the idea that 
for any truth assignment, the assignment itself or its complementary assignment satisfies at least 
one half of all clauses to multidimensional objective functions. We choose some suitable parameter 
I € N depending only on the number of objectives. For a given formula in CNF we try all possible 
partial truth assignments for each set of at most / variables using brute force and extend each 
partial assignment to a full asignment in the following way: For each remaining variable v we 
compute two vectors roughly representing the weight gained by the two possible assignments for v. 
We then compute a 2-coloring of those weight vectors with low discrepancy which completes the 
partial assignment to a truth assignment whose weight is at least one half of the total weight of the 
remaining satisfiable clauses minus some error. This error can be compensated by choosing / large 
enough such that the partial assignment already contributes a large enough weight. This results in 
a 1/2-approximation for A;-MaxSAT. 

Notations For a set of clauses H and a variable v let H[v = 1] = {C £ H \ v G C} be the set of 
clauses that are satisfied if this variable is assigned one, and analogously H[v = 0] = {C € H \v a 
C} be the set of clauses that are satisfied if this variable is assigned zero. This notation is extended 
to sets of variables V by H[V = i] = U^eV -^[^ ~ "^1 ^'-'^ i = 0, 1. 

Algorithm: Alg-k-MaxS AT (.H,w) 
Input : Formula H in CNF over the variables V = {vi, . . . ,Vm}, ^-objective weight function 

w: H ^n'' 
Output: Set of truth assignments I: y — > {0, 1} 

1 f oreach disjoint V°,V^ CV with #{V°UV^) < Ak^ do 

2 G:=H\{H[V'^ = 0]UH[V'^ = 1]); 

3 y(i-») —{veVX (yO U V^) \ 4k ■ w{G[v = i]) ^ w{H \ G)}, « = O, l; 

4 if y° n y^ = then 

5 V' ■.= v\ (yo u y^ u y° u v^), v ■.= v'yj{v\ve vy, 

6 G' := {G[V' = 0] U G[V' = 1]) \ {G[V^ = 0] U G[V'^ = 1]); 

7 for Vj G V let x^'* = E{#Ju§7) I C € G'[vj = i]} for i = 0, 1; 

8 compute some coloring x'-^' ^ {0) 1} such that 



Y^ xhxii) >- Y^ (x^'° + x^'^) - 2k6 



2 



where dr = maxjxr | Vj € V',i G {0, 1}}; 

9 let I{v) := ifoTvGV'UV'U x~^i{i}), « = 0, 1; 

10 output / 



Theorem 5.1. fc-MaxSAT is y2-approximable for any k > 1. 
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Proof. We show that this approximation is reahzed by Alg-fc-MaxSAT. First note that this algo- 
rithm runs in polynomial time since k is constant and the coloring in line [8] can be computed in 
polynomial time using Corollary 13.21 For the correctness, let (H, w) be the input where H is a 
formula over the variables V = {vi, . . . ,Vm} and w: H ^ N is the A:-objective weight function. 
Let lo'- V ^ {0)1} bs an optimal truth assignment. We show that there is a loop iteration of 
Alg-A;-MaxSAT(if , u)) that outputs a truth assignment / such that tf(/) > w{Io)/2. To this end, 
we first show that there is an iteration of the loop that uses suitable sets V^ and V^. 

Claim 5.2. There are sets V^ C I~^{{i}), i = 0,1 with #(y° U V^) < Ak^ such that for G = 
H \ (if [yo = 0] U H[V^ = 1]) and any v e V \ (F° U V^) it holds that 

w{G[v = Io{v)])<l-wiH\G). (1) 

Proof. The assertion obviously holds for ^V < 4A;^. Otherwise, we define variables {ukt+r € V \ 
r = 1,2, . . . ,k and t = 0, 1, . . . , 4A: — 1}, sequentially in a greedy fashion: 

Ho:=H 

for t:=0 to Ak-1: 
for r := 1 to k: 
j := kt + r 

choose V ^V \{ui,. . . ,Uj^i} such that Wr{Hj-i[v = Io{v)\) is maximal 
Uj := V, Hj := Hj^i \ Hj^i[v = loiv)] , Oj := w{Hj_i[v = Io{v)]) . 

Let V' = I-^{{i})n{ui, . . . ,U4fc2} for i = 0,l. We now show inequality <^, so let v € V\{V°UV'^). 
Assume that there is some r € {1, . . . , A;} such that Wr{G[v = Io{v)]) > j^Wr{H \ G). Because the 

union IJi=i ^j-A^j — ^o{uj)\ = H \G is disjoint, we get 

k 4fc-l 4fc-l 



w{H \ G) = ^ ^ akt+r' > ^ Okt+T 



r'=l t=0 t=0 

and thus 



WriG[v = loiv)]) > 



Ak-1 , . 

E \(^kt+r)r 

t=0 



Hence, by the pigeonhole principle, there must be some i G {0, 1, . . . ,4/c — 1} such that Wr{G[v = 
loiv)]) > {cikt+r)r- But sincc G C Hkt+r-i and thus even Wr{Hkt+r-i[v = loiv)]) > WriG[v = 
loiv)]) > ic(kt+r)r, the Variable v should have been chosen in step j = kt+r, which is a contradiction. 
This means that wiG[v = loiv)]) < j^wiH \ G) holds for all v G F \ (1/° U V^). D 

We choose the iteration of the algorithm where V^ and V^ equal the sets whose existence is 
guaranteed by Claim [5^ In the following, we use the variables as they are defined in the algorithm. 
Observe that by the claim it holds that loiv) = i for all v &V^ for i = 0, 1. Note that 

E ^- + ^- - E E E «^ - E #(^ n .') ji^ ^ »(G'). 

vjev vj£V'ie{o,i}CeG'[vj=i]^^ ' cgG' ^^ ' 
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Furthermore, for all Vj G V and i = 0, 1, we have the bound x-^'* < w{G'[vj = i]) < w{G[vj = 
i]) < ■^w{H \ G) because of the definition of V' and V^. By Corollary 13.21 we find a coloring 
X: y — 7- {0, 1} such that for each 1 < i < k it holds that 



2 



EE-r-E 



< 2kmax\xf''\ < 2k—Wi(H\G) = -Wi(H\G) 



and hence 

E ^''''^'^^ ^ ^ E (^''° + ^''') - ^^(^ \ G") = liHG') - w{H \ G)). 

For / being the truth assignment generated in this iteration it holds that 

wi{G € G' I I{G) = 1}) > Yj ^^'""^^'^ ^ -(^^IC) - w{H \ G)). (2) 

Furthermore, since I and Iq coincide on V \V' , we have 

w{{C £H\G' \ I{G) = 1}) = w{{C £H\G' \ Io{G) = 1}) (3) 

>w{{C(^H\G\UG) = l}) 
= w{{H\G}). (4) 

Thus we finally obtain 

w{I) = w{{C £H\G' \ I{C) = 1}) + w{{C e G' I I{G) = 1}) 

id 

> w{{G eH\G' \ I{G) = 1}) + \{w{G') - w{H \ G)) 

i w{{G eH\G'\ Io{C) = 1}) + \{w{G') - w{H \ G)) 

id 

> \w{{G eH\G' \ Io{C) = 1}) + \w{G') 
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